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The experimental availability of ultra-high-mobility samples of graphene opens the possibility to 
realize and study experimentally the “hydrodynamic” regime of the electron liquid. In this regime 
the rate of electron-electron collisions is extremely high and dominates over the electron-impurity 
and electron-phonon scattering rates, which are therefore neglected. The system is brought to a 
local quasi-equilibrium described by a set of smoothly varying (in space and time) functions, i.e. 
the density, the velocity field and the local temperature. In this paper we calculate the charge and 
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while the spin conductivity remains finite. These results pave the way to the description of charge 
transport in graphene in terms of Navier-Stokes equations. 
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I. INTRODUCTION 


Graphene, a two-dimensional (2D) layer of carbon 
atoms arranged in a honeycomb lattice, has attracted 
a huge deal of both theoretical and experimental in¬ 
terests in the last few years 1-6 . Its properties, due 
to the gapless and linear low-energy band dispersion, 
make it highly attractive for several practical applica¬ 
tions spanning the fields of optoelectronics, photonics, 
nanoplasmonics, metrology, and energy generation and 
storage 3, 1-12 . At the same time it still offers an extremely 
interesting playground for fundamental science. As the 
quality of samples continues to improve 13 , new break¬ 
throughs are expected. Many-body interactions are in¬ 
deed expected to play a crucial role in the physics of 
ultra-high-mobility samples 6 . 

The transport properties of graphene are controlled 
by the highly-mobile electrons in the n (valence) and ir* 
(conduction) bands, which arise from the hybridization of 
the p z orbitals 1 ( z denotes the direction orthogonal to the 
graphene plane). The two bands touch with a linear dis¬ 
persion at two inequivalent points ( K and K') at the cor¬ 
ners of the hexagonal Brillouin zone. For small energies, 
momenta, and doping concentrations, it is possible to ex¬ 
pand the tight-binding Hamiltonian around these points 
in a k-p fashion 1 . The resulting low-energy Hamiltonian 
describes massless Dirac fermions (MDFs) characterized 
by the density-independent Fermi velocity i>f which is 
about 300 times smaller than the speed of light. 

A high-mobility gas of free carriers can be created in 
graphene by, e.g., electrostatic gating or chemical dop¬ 
ing. We assume these carriers to be in the “Fermi- 
liquid” regime 14 . The Fermi energy is thus £f = ±vpkp 
(energies are measured from the Dirac point), where 
kp = \/ 2irn/N v is the Fermi wavevector, n is the ex¬ 
cess carrier density, and N v = 2 is the valley degeneracy. 
The sign of £f depends on the type of carriers (+ for 
electrons and — for holes). In what follows, owing to the 


particle-hole symmetry of the low-energy MDF model 1 , 
we consider only samples doped with an excess electron 
density. 

Most of the previous theoretical works on the trans¬ 
port properties of doped graphene 5 considered samples in 
which the transport is dominated by disorder effects, and 
in which electron-electron interactions play a minor or no 
role. In this paper instead, in view of the possibilities of¬ 
fered by the experimentally-available ultra-high-mobility 
samples 13 , we focus on the “intrinsic” transport regime. 
By intrinsic we mean that the transport properties are 
solely determined by electron-electron interactions, and 
disorder and phonons, as well as finite size effects, are 
considered to be irrelevant. In this Paper we focus on 
two fundamental properties of doped graphene, namely 
the charge and spin conductivities. Let us write these 
conductivities, at a finite frequency w, in the common 
form 


a^(co) 


Q 2 t Vt 

iuj + 1/tIP 


(1) 


where i = c for the charge conductivity, £ = s for the 
spin conductivity, Q c = e and Q s = h. Here r t ^ and 
are the transport relaxation times of charge and spin cur¬ 
rents, respectively, while V c and V s are the corresponding 
“Drude weights”. 

Any conductivity, associated with the transport of a 
physical quantity, are, in general, affected by electron- 
electron interactions. Both its “Drude weight” and its 
relaxation time change when interactions are turned on. 
In fact, as it happens for the quasiparticle lifetime (r ee ), 
also transport relaxation times are limited by thermally- 
activated electron-electron scattering processes. More¬ 
over, interactions renormalize the quasiparticle proper¬ 
ties as the Fermi velocity, and accordingly the Drude 
weights. The self-energy corrections which “dress” the 
interacting Fermi velocity, Vp, are especially important 
in graphene, since its low-energy MDF Hamiltonian has 
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an infinite bandwidth 1,6 . Interactions between electrons 
at the Fermi surface and states at large negative ener¬ 
gies (in the valence band) lead to a logarithmic diver¬ 
gence of 4 15 as the system approaches the undoped 
regime. Furthermore, electron-electron interactions are 
also responsible for “vertex corrections” to the Drude 
weights, which are encoded in the Landau parameters 14 
Fn s . Vertex corrections are usually small in a wide 
range of interaction strengths 14 . The full renormaliza¬ 
tion of the charge Drude weight takes the form 14,16,17 
V c = T>i°^Vp(l + F{)/v f, while the spin Drude weight 
becomes 18 V s = + F*)/v f (see also App. F). 

Although this is a very general scheme, care must be 
exerted when dealing with charge currents. In the case 
of a parabolic band 2D electron gas (2DEG), for exam¬ 
ple, Galilean invariance leads to a perfect cancellation 
between self-energy and vertex corrections to the charge 
Drude weight. In formulas ^(l+F-j 5 ) = Vy. For the same 
reason, the charge transport time of a clean Galilean- 
invariant Fermi liquid is unaffected by electron-electron 
interactions in the absence of Umklapp processes, and 
it is thus infinite as in the non-interacting case. These 
facts show that electron-electron interactions, which con¬ 
serve the total momentum of the system at any scatter¬ 
ing event, are also inefficient in relaxing a homogeneous 
current. In Galilean invariant systems, the latter is in¬ 
deed proportional to the total momentum, which is a 
conserved quantity. At odds with this, the spin current 
follows the “general rule” outlined above and is relaxed 
by electron-electron interactions, which can transfer mo¬ 
mentum between the two spin populations, giving rise to 
the phenomenon of the spin-Coulomb drag 19-21 . 

Since graphene is not a Galilean invariant system, 
charge relaxation due to electron-electron interactions is 
not forbidden a priori by any symmetry. In the low- 
energy MDF model, indeed, the current and total mo¬ 
mentum are not proportional to each other. It is there¬ 
fore not surprising that electron-electron interactions af¬ 
fect the charge Drude weight of graphene in a non-trivial 
way 16,17 , i.e. the product Vp(l + Ff) ^ vf- It has been 
shown that, to the first order in the strength of electron- 
electron interactions 16 , vertex corrections exactly cancel 
the self-energy renormalization due to particle-particle 
scattering at the Fermi surface, but do not affect the log¬ 
arithmic divergence of the Fermi velocity 16 . Moreover, 
while it is clear that in a 2DEG the charge transport 
time is unaffected by interactions, no conclusion can be 
drawn a priori for graphene. 

In this paper we prove that at low temperature the 
charge transport time is infinite, i.e. 1 /r^f 1 = 0 + , while 
the spin transport time r t ; r 1 is limited by electron-electron 
interaction and it is thus finite. As noted above, the rela¬ 
tion between the current and total momentum is highly 
non-linear in graphene. Is it thus somewhat surprising 
that the charge transport time is not affected by electron- 
electron interactions, and that only the Drude weight is 
renormalized. This result can be understood as follows. 


While the momentum k and velocity v\ = XvFk/\k\ 
of a quasiparticle are not directly proportional to each 
other, they become approximately linearly related at low 
temperature for any finite doping concentration. Indeed, 
in the limit of k^T -C £f, the dominant contribution 
to the transport comes from electrons in a thin shell of 
size k^T around the Fermi energy. All these quasipar¬ 
ticles have magnitude of momentum equal to kp, and 
velocity 1 ?+ ~ vpk/kp , if the system is n-doped. This 
in turn implies that a linear relation is established be¬ 
tween the momentum and velocity of each quasiparticle 
and, accordingly, between the current and total momen¬ 
tum of the system. Since the latter is conserved, at low 
temperature electron-electron interactions cannot relax a 
homogeneous current and doped graphene behaves as an 
effectively Galilean-invariant system. 

We stress that this argument applies only to the calcu¬ 
lation of the charge transport time and breaks down when 
one considers the charge Drude weight. The latter has 
contributions from virtual processes between all quasi¬ 
particle states, not only those around the Fermi energy. 
Since these processes span all the quasiparticle spectrum, 
the non-linear relation between the current and momen¬ 
tum operators becomes apparent and the Drude weight 
gets renormalized. 

The situation is completely different for the spin con¬ 
ductivity. In this case electron-electron interactions (i) 
renormalize the spin Drude weight and (ii) provide a fi¬ 
nite transport time for spin currents. While the former 
effect is expected to be small in a wide range of values of 
the strength of electron-electron interactions 18 , the lat¬ 
ter is large. The spin conductivity, which was infinite 
in the non-interacting limit, turns out to be finite in 
an interacting system. It thus offers a more powerful 
probe of electron-electron interactions as compared with 
the charge conductivity. The physics behind the spin 
conductivity is intimately related to the phenomenon of 
spin-Coulomb drag. When a pure spin-polarized current 
is injected into the system, each spin component of the 
current exerts friction on the other spin component via 
electron-electron interactions. The relative velocity of 
the spin populations therefore decays in time, and even¬ 
tually vanishes unless an external driving field is present, 
in which case it reaches a steady state. This in turn im¬ 
plies that the spin conductivity must be finite. Our cal¬ 
culation shows that typical values for the spin transport 
time range between 1 — 10 ps. 

This paper is organized as follows. In Sect. II we de¬ 
fine the low-energy MDF model of graphene, and we set 
up the all-order diagrammatic calculations needed to de¬ 
termine the charge and spin conductivities. The main 
steps of the calculation are given in Sect. Ill, which also 
presents the main results of our paper, namely the charge 
and spin transport times. Our results are summarized in 
Sect. IV. Appendices. A-F provide several technical de¬ 
tails of the calculation. 
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II. MODEL AND BASIC DEFINITIONS 


We model graphene with the low-energy MDF Hamil¬ 
tonian (per valley flavor - hereafter h = l) 1,6 

Ti = Y £ k,\^Lx,a^k,\,a + ^Y W*-<j ~ ™o) , (2) 

k, X q 


where ipk,x,a (V'l a <t) destroys (creates) a particle with 
momentum k and spin a = ± in band A = ±, £k,x = 
Xvpk, and v q = 2ne 2 /(eq) is the non-relativistic Coulomb 
interaction. Here e models the dielectric environment 
surrounding graphene and, as a first approximation, it 
is the average of the dielectric constants of media above 
(ei) and below (€ 2 ) the sheet, i.e. e = (ei + £ 2 )/%. The 
strength of electron-electron interactions is characterized 
by the density-independent “fine-structure constant” of 
graphene (restoring h) a ee = e 2 / (hev f). Finally, the den¬ 
sity operator is 

fl q= Y v ^'( k -’ k +)^k_,x .o-V’fc+AV , (3) 

k,o, A,A' 


where we defined k± = k ± q/2 , and the matrix element 
of the density operator between the eigenstates of the 
bare Hamiltonian is 1 

i(<p k -ip k ,)/2 , \\' e -i((pk-Vk')/2 

Vxx<{k,k') = - - - . (4) 

Here < pk is the angle between the momentum k and the 
cc-axis. 

The spin-resolved current operator of MDFs is 

3q ] = Y J ^y( k -> k +)^k_,X,a^k + ,X’,a , ( 5 ) 

k,\,\' 

where Jxx'(k,k') = [J^, (fc, fc'), (fc, fe')], and 

, , Vpi-iVk+Vk ’)/ 2 \p-i(Vk+Vk’)/ 2 

- , 

, , VpiiVk+Vk')/ 2 — \ e -i(Vk+Vk’)/ 2 

J AA'( fc , fc ') = V F - ’( 6 ) 


are the matrix elements of the current operator between 
the eigenstates \k, A) and |fc', A') of the system. 

Fig. 1 summarizes the all-order diagrammatic re¬ 
summation needed to calculate the charge and spin con¬ 
ductivities, 22,23 which are microscopically defined as 24 




26 v—t . x 

OJ z — J a Jp 


( 7 ) 


and 


(s) 
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z j p 

cro' 


( 8 ) 



FIG. 1: a) The diagrammatic representation of the current- 
current response function. The left dot is the bare vertex 
A (°’“) (we suppress the momentum-energy dependence for 
brevity), while the solid double lines are Green’s functions 
dressed by the self-energy. In the large-np limit it correspond 
to the GW self-energy, which is depicted in panel b). Wavy 
lines represent RPA screened interactions. Finally, the tri¬ 
angle represents the vertex function A' 3 which is dressed by 
e-e interactions and satisfies the Bethe-Salpeter equation in 
panel c). Note that the form of the irreducible interaction I is 
uniquely determined by the choice of the self-energy, provided 
A' 9 must satisfy the Ward identities 14 (see Fig. 2). 



FIG. 2: The diagrams that contribute to the irreducible 
interaction I of Fig. 1. 


Here X .(„) .<.*’) {q = 0,w) is the proper spin-resolved 

Jp 

current-current linear response function 14 given by the 
diagram in Fig. la). Its analytical expression on the 
imaginary-frequency axis reads 

Kj (?) *^m) = N v kftT ^ ^ i ''a (k—,i£ n ) 

k,X,X' £ n 

x A^r ) Q (fc_,fc + )G^ ) (fe+,i£n+*w m ) 

x ^yx,l( k +’ i£n + iu; mi k— i i£ n ) ■ (9) 

Here a, (3 = x,y denote the Cartesian components of 
the vectors. The double lines in Fig. 1 correspond to 
the spin-resolved Green’s functions which are dressed by 
the GW self-energy insertion of Fig. lb) (see Sect. IIA 
below). In Eq. (9) £ n = (2 n + 1 )/fi (oj m = 2m/f3) 
are fermionic (bosonic) Matsubara frequencies (n, m = 
0,1,...), G^\k,i£ n ) is the Green’s function on the 
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imaginary-frequency axis, while the bare current vertex 
determined from Eq. (5) reads 

A^ a (k-,k + ) = j[ a J(k_,k + ). (10) 

The term in the last line of Eq. (9), namely 

A*va \(k+,ie n + iui m ,k-,i£ n ), is the vertex function, 
which is dressed by electron-electron interactions and 
satisfies the self-consistent Bethe-Salpeter equation of 
Fig. lc) (see Sect. IIB below). We stress that the choice 
of the self-energy, together with the requirement of ful¬ 
filling the Ward identities, uniquely determines the self- 
consistent Bethe-Salpeter equation satisfied by the vertex 
function, i.e. the irreducible interaction I. 14 
The GW self-energy shown in Fig. lb) reads 

E ^\k,ie n ) = -fc B T ^ ^ W(k' - fc, is n < - ie n ) 

k', X' e n f 

x G { y{k',ie nl )V xx ,{k,k')V x , x (k',k) . 

( 11 ) 

Here is the screened electron-electron inter¬ 

action, represented in Fig. lb) by a wavy line. In the 
large-IV v limit this is given by 


The three contributions Ay’^ a /3 ' 1 (k + ,ie n + fc_, «£„) 
(with i = 1,... ,3) correspond to the three diagrams in 
Fig. 2. They read 


Aftjr ) ( fe +> + iu} m, k- , i£„) = - k B T ^2 

k',e n / /ij/i/jCr" 

xW \\\k',k,i£ n i - ie„)G^, \k' + ,i£ n : + iu] m ) 

xG^'^(k'_,i£ n >)A { 2{k' + ,i£ n ' + iu m , k'_, ie n >) , 

(15) 


and 


■^A'A ,(3 (^+ ? ^mi 5 ^n) — ^ ^ 

k' ,En' 

x W xxw ‘ ) ( fc,; fc ’ i£n ' + i£n + iuj ™) 
xG { 2 \k' + ,is n > +iu m )G ( £ \k'_,i£ n f) 

xA ^ ) ( A;, + ’* e "' +iu} mi k— 1 l£ n ') ■ (16) 

Here we define 


W(q,ifl m ) 


V q 

1 ^ qXnniQ 7 ^m) 


( 12 ) 


where Xnn(q,u) is the proper density-density response 
function 14 of graphene. In principle, this should be cal¬ 
culated in analogy to the current-current response func¬ 
tion of Eq. (9), i.e. it should contain dressed Green’s 
functions and vertex function. However, to simplify our 
calculation we neglect the vertex corrections to Xnn(q , u>), 
which we define as 

Xnn{q,iu m ) = N v k B T ^ ^ G^n } (<?', lEn) 

q' ,En A " ,1-L" 

x G^, } (q' + q,i£ n + iui m ) 
x £>a>" (<?', q' + <?)1V'A" (<?' + <?,<?') • 

(13) 


W \x2n'\ k '^ k '< iuJ rn) = 5 a< j"W(k - k', iuj m ) 
xD AV (fc + ,fc' + )^(fe'_,fc_) , (17) 

k ’ i£ n' - i£ n) = N v k B T J2 

x W(q',iuJm’)W(q'- q.iUm’~ iu m ) 

\",n" 

xVy X "{k + ,k + - q')V X n X (k+ - q ', fe_) 

xT> w //(fc , _,fc' + - <z')2Vv( fe + - 

xG a „ (A;_[_ q , T 2cc m iuj m '^ 

xG^„ (&_{_ q :i^n' T i^m ') j (13) 


Note that the density vertices in Eq. (13) are not renor¬ 
malized by electron-electron interactions. This approx¬ 
imation does not invalidate the results obtained from 
Eq. (9), i.e. it does not lead to qualitatively wrong behav¬ 
iors of the conductivities. This because both self-energy 
and vertex corrections are taken into account on an equal 
footing in Eq. (9). The approximation done in Eq. (13) 
has however an impact on the form of the Bethe-Salpeter 
equation satisfied by Ay X p(k + ,i£ n + iuj mi fc_, i£ n ), i.e. 
on the form of the irreducible interaction I [see Fig. 2], 
Indeed, the requirement of fulfilling the Ward identities 14 
constrains the dressed vertex to satisfy the following self- 
consistent Bethe-Salpeter equation 14 : 

Aa'A,/?(^+> n ^mi k —, i£ n ) A A /a(&+, k— ) 

+ A x ; x a p\k + ,i£ n +iu) m ,k~,i£ n ) . (14) 

2 = 1 ,.. .,3 


and finally 

W >of2' ) ( k '’ k ’ i£ n' + ie n + iuJ rn ) = N v k B T ^ 

q' ,cJ m / 

x ^ W(q',iu m ')W{q' -q,iu m' ^m) 

X" 

xV xx »(fe_,fc_ + q')V X n X ,{k- +q',k +) 
xV^»(k'_,k' + - q')£Vv( fe + “ q\k' + ) 
x Gy) (fe_ + q', i£ n + i<jj m >) 
xG- q',i£ n ' +*w m - iu m >) . (19) 

In what follows we start from the evaluation of the self¬ 
energy corrections, and of the quasiparticle lifetime at the 
Fermi surface, and we then proceed to the calculation of 
the vertex correction. 
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A. The quasiparticle decay rate 


In this section we calculate the quasiparticle lifetime 
r ee at the Fermi energy, defined as 


1 

Ae 


OS 


V CT, 0M + )] 


k—kr. 


,( 20 ) 


where the self-energy was defined in Eq. (11). We focus 
on the imaginary part of the self-energy, which controls 
the charge and spin relaxation times, and we disregard 
its real part, which is responsible for the renormalization 
of the Drude weights. In the spirit of Landau theory of 
normal Fermi liquid, we take care of this approximation 
of the diagrammatic calculation by replacing a posteri¬ 
ori the non-interacting Drude weights with their inter¬ 
acting values. A microscopic calculation of the charge 
Drude weight to the first order in the strength of electron- 
electron interactions was given in Refs. 16. As shown in 
App. A, at low temperature (ksT <C £f) Eq. (20) be¬ 
comes 


qp 


(fc B T) 

37r 


2 r 2k F (l-T/2) 


'k F T 


dq q\W {q 7 Tep)\ 2 


X A(,,-f ep) , 


( 21 ) 


where T = fkp,T/e p, £ = 7r/-\/5, and 


7T f 2 * 

A(q,w) = -/ dip„5(\k — q\ — k — uj/vp) 

«F J o 


«F Jo 

1 + COS {iPk-q + Vk) 


477 


k—kp 

kp — Lj/vp 


V F ^/^k^q 2 — ( q 2 — l o 1 1v\ + 2kpco/vp) 2 

1 q^^/vl \ 


x 0 1 - 


4kp(kp — u/vp) J 

q 2 — uj 2 /vp + 2kpto/vp 


2kpq 


( 22 ) 


Eq. (21) describes, as shown in Fig. 3, the decay (scat¬ 
tering) of a quasiparticle of momentum k to a state of 
momentum k — q through the creation of an electron- 
hole pair of total momentum q obtained by exciting a 
particle of momentum k" — q to a state of momentum 
k". Such a process is encoded in the density-density re¬ 
sponse function 25-29 STOx nn (q, Tep) and is depicted in 
Fig. 3. Notice that, since all the initial and final states 
are on the Fermi surface, the conservation of momentum 
implies that k and k" — q (and thus k — q and k") are 
diametrically opposite. This fact will be used in what 
follows to simplify the expressions of the transport times 
(see App E). 

Numerical results obtained from Eq. (21) are shown 
in Fig. 4. In passing, we recall that the GW quasiparti¬ 
cle lifetime is also calculated in Refs. 30,31. Contrary to 



FIG. 3: A pictorial representation of double particle-hole 
excitations that contribute, to lowest order in the strength 
of e-e interactions, to the quasiparticle decay rate calculated 
in Sect. II A. Note that, since all the states involved in the 
scattering process live at the Fermi surface, the conservation 
of momentum constrains the initial states k and k" — q to be 
diametrically opposed. The same happens to the final states 
k — q and k". 


Eq. (21), the expressions given in Ref. 30 do not contain 
any low-temperature approximation (see App. A for more 
details). In Figs. 5-6 we show a comparison between the 
quasiparticle lifetime calculated from Eq. (21) and the 
“exact” one computed in Ref. 30. Note that the agree¬ 
ment is very good in the chosen range of temperatures 
and densities. 


B. The vertex function and the Bethe-Salpeter 
equation 

In this section we summarize the complicated calcula¬ 
tion of the vertex correction to the charge and spin con¬ 
ductivities. The details of the derivation can be found in 
Apps. B and C. 

The first step of the calculation is to analytically con¬ 
tinue Eq. (9) to real frequencies. We indeed recall that 
the charge and spin conductivities are determined, ac¬ 
cording to Eqs. (7) and ( 8 ), by the u —> 0 limit of the 
retarded spin-resolved current-current response function 
calculated at q = 0 . We stress that the analytical contin¬ 
uation ico m —» uj + ir\ (with q = 0 + ) must be performed 
before the small-frequency limit. 

After the analytical continuation to real frequen¬ 
cies 32 , the current-current response function of Eq. (9) 
contains products of advanced-advanced (schemati¬ 
cally G a G a ), retarded-retarded (G R G R ) and advanced- 
retarded (G a G r ) Green’s function. The first two contri¬ 
butions, G a G a and G R G R , have poles on the same side 
of the complex plane. In the limit £FT ee 1 we can ne¬ 
glect them 32 and retain only the “mixed” term G A G R . 
After some simple algebra, shown in detail in App. B, 















6 




T[K\ 


FIG. 4: (Color online) Panel a) the quasiparticle lifetime of 
massless Dirac fermions r eB , as defined in Eq. (21), in units of 
picoseconds and plotted as a function of the density n in units 
of 10 12 cm -2 for three values of the dimensionless coupling 
constant a ee . In this plot we fixed the temperature T = 
300 K. Panel b) same as in panel a) but shown as a function 
of temperature (in units of K) for a fixed excess carrier density 
n = 10 12 cm -2 . 
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Eq. (9) becomes 

/ de 

— [n F (e + w) — n F (e)] 

K,A,A- m 

x G[ A ’ a) (fc_, e)A, K)G { * y) (k + ,6 + w) 
x A va^( fe+,e++w,fc-,e-) . (23) 


FIG. 5: A comparison between the quasiparticle lifetime 
calculated from Eq. (21) and the one computed in Ref. 30. In 
this figure the temperature is kept fixed at T = 300 K and the 
quasiparticle lifetime is plotted in units of ps as a function of 
the density (in units of 10 12 cm -2 ). Panel a)-c) refer to the 
three values of the coupling constant a ee = 0.5, a ee = 0.9, 
and a ee = 2.2, respectively. 


1,..., 3} on the last line of Eq. (24) read 


Accordingly, after the continuation to real frequencies 
the self-consistent Bethe-Salpeter Eq. (14) becomes 


A VAj( fe +> e + +w,fc—,£_) = <W A^^(fe+,fe_) 

+ A a'’a,/3 \k+,e + + u>,k_,E-) . (24) 


Also in the derivation of Eq. (24) we retained only mixed 
terms of the form G A G R . In spite of this approxima¬ 
tion, our calculation yields a closed set of self-consistent 
equations. The terms { A vaV (^+> £ + = 


JV A'A,/3 


{k+,£+ + w,fc_,£_) = 47V v ^ ^2 ^2 


k' ,q' ,a" X" ,/x" 


x [n F (e') + n B (£' - e)] [n F (o/ + s') - n F (cj' + e)] 


xSm 

xSro 


G^y,’ a \q' — k,uj' + e) 

G { y\q'~k'^+e') 


'X-'Dx'fi' {k .|_, k + )D ll \(k _, fc_) 

xT>a'>"(<?' ~k,q' - k') r D ll >'x"{q' - k',q' - k) 

xG^’^(fc' + ,E' + ui)G^ ,rT \k'_,e') 

x A ^>i( fc +> e + + fe -’ e -) > 


(25) 
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and finally 

A v^ ) ( fe +' £ ++ w ' fe -- £ -) = - 4Ar v 

k' ,q' ,<j" X" 

x /£/S 

x [n F (e') + n B (£' + e)] [n F (u/ + e) — n F (o/ - e' - w)] 
xSm[G^ ,,T '(fc_ + qr', e + u/)] 
x3?n[G^’ a \k' + — q', e' + ui — <*/)] 
xD AA //(fc_, fc_ + q')V\n X ,(k_ +q\k+) 
'x'D^”{k'_,k' + - q')V^^f(k' + - q\k' + ) 

xG^"\k' + ,e' + uj)G^"\k'_,e') 

xA >>,V(^+’ e + + w ’ e -) ■ (27) 

Equations (23)-(27), together with the inverse quasipar¬ 
ticle lifetime defined in Eq. (21), constitute a closed 
set of equations that can be used to determine the re¬ 
tarded spin-resolved current-current response function in 
the limit v F q <C uj,t~,} <C £f- This calculation is per¬ 
formed in the next section. 


III. THE CHARGE AND SPIN 
CONDUCTIVITIES OF GRAPHENE IN THE 
FERMI-LIQUID REGIME 


FIG. 6: A comparison between the quasiparticle lifetime 
calculated from Eq. (21) and the one computed in Ref. 30. In 
this figure the density is kept fixed at n = 10 12 cm -2 and the 
quasiparticle lifetime is plotted in units of ps as a function 
of the temperature (measured in K). Panel a)-c) refer to the 
three values of the coupling constant a ee = 0.5, a ee = 0.9, 
and a ee = 2.2, respectively. 


and 


In this section we present the derivation of the charge 
and spin conductivities of graphene. We start from 
Eqs. (23)-(27), where we set q = 0 and we take the limit 
w —» 0. To O(uj) Eq. (23) becomes 


X,W,(. r'){q = 0,w) = U}N V V 

Joe Jp Z ' 

fc, A,A' 


de / dn F {£)\ 
27ri \ de ) 


x G{ A ’ a) (, k , e) A ( ^ a (k, k)G^’ a) (k,e + uj) 
x A< x'x}(k,e + +uj,k,e-) . (28) 


a v A ™0 ) (k + ,e + +u>,k-,e-)=4:N v H H 

k' ,q' ,cr" “ >>" 

x [n-p(e') + 71b(£ / - £)] 
x \n F (uj' — £ — uj) — rip (a/ — s’ — w)] 
x$sm[G^,’ <J \k + — q 1 ,e + uj — w')] 
x3m[G^’ cr \k' + — q' ,e'+ uj — uj')] 
xT> X ’\"(k + ,k + - q')V X "x(k + - q',fc_) 
xV^" ( k '_, k’ + - q')l(k' + - q\ k' + ) 

xG^"\k' + ,e' + u,)Gl A ’°"\k'_,e') 
xA ^ M ,V( fe +> £ + + w,fc'_,e'_) , 


The limit w —> 0 is understood in this equation. We 
observe that the function —dn F {e)/de is peaked around 
£ = 0 and tends to a 5-function in the low-temperature 
limit £ F r ee 1. We thus evaluate all the other functions 
on the right-hand side of Eq. (28) at £ = 0, with the 
exception of A^7 a ^(fe, £+, k, £_), which requires further 
care. As it will become clear in what follows, the latter 
contains Fermi and Bose factors which depend on £ and 
that combine with dn F (e)/de in Eq. (28) to yield the 
correct transport times. Missing this step would lead 
to a non-cancellation between the self-energy and vertex 
corrections in the charge channel. In the limit w-)0we 
also approximate 


g(A,dr) 


(fc,0)Gi R ’ CT) 


(fc,w) 


2 iSxx' 
uj + i/r ee 


[Gf’ CT) 


(M)] • 

(29) 


(26) 













In so doing we neglect the incoherent part of the 
Green’s function, i.e., the part of G that is not in¬ 
cluded in the quasiparticle-pole approximation. Herein 
lies our Fermi liquid approximation. At low temper¬ 
ature, 3m[G^ R) (fc, 0)] is a Lorentzian strongly peaked 
around £k,\ = £f and with a width proportional to the 
quasiparticle decay rate. This implies that k ~ ky and 
A = A' = +. For £FT ee 1 the transport is dominated by 
states that lie in a thin shell of thickness ~ k B T around 
the Fermi energy 14 . 

At the same level of approximation, the dressed vertex 
function satisfies the following Bethe-Salpeter equation 
(see also App. D): 


A++^(fe,£+ + w,fc,e_) = <WA++°^(fc,fc) - 


(0,cr) 


8 iN v 


+ i/i 


k' ,q' ,a" 


x \ny{uj’ + s') — ny(uj' + £)]3m[G^ R,<T \k'. 0)] 
x 3m [G^ 1 ’^ (fc — q ', 0)]Sm[G! | R ' ff \q' — k! , 0)] 
xV ++ (k, k - q')V ++ (k - q ', k)V ++ {k ', fc' - q') 


xV ++ (k' -q',k') 


A 


(aa) 


(fc - q , e, + w, k - q ', e '_) 


+A ^\k\e' + +u,k',e'_) 


++,/3 \ K 

-A + + a B \k' - q',e' + +uj,k' - q',e'_) 


(30) 


Here the limits k = ky and u> —> 0 are understood. In 
deriving Eq. (30) we used that in the low-temperature 
limit also the momenta fc', k — q' and k' — q' are all 
pinned at the Fermi surface, i.e. |fc'| = \k — q'\ = 
\k' — q'\ = ky. and that the corresponding quasiparti¬ 
cles live in the conduction band (recall that the system 
is n-doped). Thus Eq. (30) is a closed self-consistent 
equation for the dressed vertex with all momentum- and 
energy-arguments pinned at the Fermi surface. 

Equation (30) can be further simplified by noting that 
its solution must be introduced into Eq. (28). We can 
thus immediately carry out the angular integration over 
£ with the weighting function dny(£)/d£, and we get 


A++J,(fc,u; + ,fc, 


|W(^0)| 2 3m[Gf’°(fe',0)] 

i [GV R ’ CT) (fe - q', 0)] 3m [G^° {q' - 
+ (fe, fe ~ < 7 , )'C , ++(fe - q', k)V ++ (k' 

/» / / * t\ i (era') / * / _L- i 


x3?n 
xV ++ (k 
xP + ^fc'-q',fc')[A^}(fc-q', 

+A^' ) (fc , , w +,fc , ,0-) 

-A+ + <T « ) (fc' - - q',0~) 


A ++}(k - q'.W 


-A + + a ft\k' - q',w + ,i 
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To get this equation we used the results of Eqs. (A12)- 
(A14). 

We now solve Eq. (31) with standard methods 32 . We 
first reduce it to an algebraic equation with the following 
Ansatz: 

A++J,(fe,w + ,fc,(T) = 7 o-<t' (v)A++J (fe, k) , (32) 


and we then solve it for 7 CT<T '(w). Note that at the Fermi 
surface A ^ 0 )) 7 2 (fc, fc) = vykp. A further simplification 
comes from the fact that the Green’s functions on the 
right-hand side of Eq. (31) are independent of spin. We 
suppress their spin dependence in what follows. 

Let us first consider the Bethe-Salpeter Eq. (31) in 
the charge channel. To obtain the dressed charge vertex, 
after having introduced the Ansatz (32) into Eq. (31), we 
sum over the spin index a. After some algebra we obtain 


[l+a' M + 7-rr' M] A^j (fe, fe) = A ( °fJ (fe, fe) 

~ 3("+Vr.|) [7+J ' M+7 ^' Ml £ IWW’ 0 )! 2 


x3m [Gf (fe', 0)] 3m [G^ R) (fe - q', 0)] 
x3m[Gf(g' - fc',0)]D ++ (fc, fc - q') 
xV ++ (k - q', fc)2> ++ (fe', fe' - q')v ++ (k' - q', fc') 


A^+J(fc - q', fe - q') + A^J(fe', fc') 
-A ++}{k' - q', fc' - q') 


(33) 


Let us now consider Eq. (31) in the spin channel. In this 
case we first multiply by a and then we sum over the spin 
index a. We get 


\l+a' M - 7-ct' M] A^j (fc, fc) = a' (fc, fc) 
8iN v (k B T) 2 f ^ ^ | w, i m |2 

“3(J+i7^) [7+ "' ( )] ^ 1 ta ' 0)1 

v ' k' ,q' 

x3m[G™(fc',0)]3m[Gf h R) (fc - q',0)] 
x3m[G™(g' - fc',0)]2?++(fc, fc - q') 
xV ++ (k - q', fc)27 ++ (fc', fc' - q')Z>++(fc' - q', fc') 


/ A I 0 ’' 7 ') 


xA^j(fe - q',fc - q') . 


(34) 


Eq. (33) and (34) control the vertex renormalization in 
the charge and spin channels, respectively. Note that the 
former depends only on the variable 7 +£r /(u;) + 7 _ CT '(a;), 
while the latter depends only on 7 +( 7 /(u;) — 7 -^/(a;), and 
that they are thus independent of each other. The spin 
and charge channels are thus decoupled in an unpolarized 
system. 

We are now in the position to solve Eqs. (33)-(34). We 
project them along the direction of A^fi (fc, fc) = vykp 
and, after some straightforward algebraic manipulations, 
we find 


7+<t'M +7-<t'M 


W -(- i / T ee 


(31) 


OJ + IT] 


(35) 













9 




T[K\ 


FIG. 7: (Color online) Panel a) the spin transport time 
of massless Dirac fermions as defined in Eq. (37), in 
units of picoseconds and plotted as a function of the density 
n in units of 10 12 cm~ 2 for three values of the dimensionless 
coupling constant a ee - In this plot we fixed the temperature 
T = 300 K. Panel b) same as in panel a) but shown as 
a function of temperature (in units of K) for a fixed excess 
carrier density n = 10 12 cm -2 . 


and 


7+cr' (w) - 7-0-' M = <?' U + l ^T s) ■ (36) 

We recall that the quasiparticle lifetime at the Fermi sur¬ 
face, r ee , is defined in Eq. (21) and it is explicitly calcu¬ 
lated in Fig. 4. In Eq. (36) we defined 

4 ) = -\K{k B TfY J \ W ^'M 2 ^rn[Gf\k-q',())\ 

^tr 

x ^sm[Gf-\k',0)]^sm[Gf\q' - fc',0)] 
x X> ++ (fe,fc- q')T> ++ (k- q ', fe) 
x V ++ (k',k'- q')V ++ (k'- q',k') 
x [l - cos(</? fc - Vk- q >)\ ■ (37) 

The calculation of Eq. (37) is performed in App. E. Nu- 

(s) 

merical results for rf r are shown in Fig. 7. 

We substitute Eqs. (32), with 7 CTCT /(u;) defined in 
Eqs. (35) and (36), back into the definition of the spin- 
resolved current-current response function of Eq. (28). 
From this we then compute the charge and spin con¬ 
ductivities, according to the definitions given in Eqs. (7) 
and (8). After some lengthy but straightforward algebra 


we get (we restore H in the following expressions) 


r (c) ( 
7*7 \ 




X 

X 


n -t £ £ - 


h ' / 27T 
fc.A.A' ’ 

A°) (U 


de f dn F (e) 
de 


G{ A) (k,e) 


AZ'Jk, k)GZ (fe, £ + w)A<J? Ai/J (fe, fe) 

'W'M 


a,a '=± 
&a{3 


e 2 V { Z 

—iuj + 77 


(38) 


In Eq. (38) we defined the non-interacting charge Drude 
weight T>Z = N v £F/(2nh 2 ). Note that in this equation 
the quasiparticle lifetime disappears as a consequence of 
a cancellation that occurs between the product of the two 
Green’s functions, approximated as in Eq. (29), and the 
vertex correction (35). 

Equation (38) shows that the real part of the DC 
charge conductivity of interacting graphene is infinite. 
Moreover, the weight of the low-frequency Drude peak 
coincides with the non-interacting one. As discussed in 
Sect. II A, our theory does not capture the renormaliza¬ 
tion of the Drude weight due to electron-electron inter- 
actions 6 ’ 16,17 , since we have neglected the contribution of 
the real part of the one-body self-energy. Such contribu¬ 
tion can be taken into account by replacing 


Vi 0) -¥ V c = 


W 

2irh 2 


(1 + F*)vpk F . 


(39) 


Note that electron-electron interactions affect the Drude 
weight in a two-fold way. On one hand the bare Fermi ve¬ 
locity is replaced by its interacting counterpart Up, which 
embeds the self-energy corrections at the single-particle 
level. On the other hand, electron-electron interactions 
provide also “vertex corrections” to two-body properties. 
In the case of the charge conductivity these are encoded 
in the Landau parameter 14 F®. Both the self-energy and 
vertex corrections to the charge Drude weight have been 
calculated to the first order in the strength of electron- 
electron interactions in Ref. 16. 

Our calculations show that doped graphene at low tem¬ 
perature behaves as an effectively Galilean-invariant sys¬ 
tem. Indeed, if graphene is doped and the tempera¬ 
ture is sufficiently low, the velocity and momentum of 
the current-carrying states, i.e. the quasiparticles at the 
Fermi surface, are linearly related via the Fermi wavevec- 
tor fep- This in turn implies that the current carried by 
such states is proportional to their total momentum. As 
in a parabolic-band electron gas, electron-electron inter¬ 
actions conserve the total momentum and are thus in¬ 
efficient in relaxing the current carried by quasiparticles 
at the Fermi surface. We emphasize that this conclusion 
is reached for doped graphene. As the doping level de¬ 
creases at a given temperature the momentum-velocity 
relation can no longer be linearized and the behavior of 
the conductivities must be reconsidered. This is in all 
likelihood the reason why pristine (undoped) graphene 
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exhibits a finite charge conductivity, as opposed to the 
infinite conductivity discussed here. 

Following steps analogous to those taken in the deriva¬ 
tion of Eq. (38), we get the spin conductivity 


(s) 

T ap 


M 


X 

X 



A^\k,k)G^(k,e)A^f(k,k) 


aa ' ’Ter cr' (fcp) 

cr,(T / =± 


^Ck/3 


h 2 vi 0) 

-iu + l/r t ( r s) 


(40) 


Where the non-interacting spin-Drude weight vi 0 ' 1 coin¬ 
cides with T>i°\ Again, our result misses the renormal¬ 
ization of the spin-Drude weight, which is calculated in 
Ref. 18 for a parabolic-band electron gas and in App. F 
for graphene and reads 


V s 


1 + Ff 
1 + Ff 


V c . 


(41) 


As it happens in a 2DEG, this renormalization is ex¬ 
pected to be small also in graphene, for typical carrier 
densities and in a broad range of values of the coupling 
constant a ee . The finiteness of the relaxation rate is thus 
by far the largest effect of electron-electron interactions 
on the spin transport in graphene. The spin conduc¬ 
tivity becomes finite when particle-particle interactions 
are turned on, even in a perfectly clean and infinite sys¬ 
tem, and scales as ~ T~ 2 in the DC limit (to —> 0). 
Finally, we can calculate the spin-drag transresistivity 21 
Pl! in the limit kpT -C £p, according to the defini¬ 
tion p-jq = — [r t ( r s) T> s /e 2 ] -1 . Neglecting the Fermi-liquid 
renormalization of the Drude weight, For k^T/ep = 0.1 
and for a coupling constant a ee = 2.2 we get |p^| ~ 54 D. 
As a comparison 21 , in a 2DEG with r s = a ee and same 
ratio kpT/ep, Ip-jql ~ 200 D. 


IV. SUMMARY AND CONCLUSIONS 

In this Paper we have calculated two fundamental 
transport properties of doped graphene, the charge and 
spin conductivities, in the “hydrodynamic” regime of 
the electron liquid 33 — see Eqs. (38) and (40). In this 
regime the rate of electron-electron collisions is extremely 
high, and the other scattering mechanisms, like electron- 
phonon or electron-impurity interactions, are negligi¬ 
bly, 35 Electron-electron collisions bring the system to 
a state of local quasi-equilibrium, which is described by 
a set of smoothly varying (in space and time) functions, 
i.e. the density, the velocity field and the local tem¬ 
perature 36 . Such intriguing regime is rarely relevant in 
solids 33 , since momentum-non-conserving processes dom¬ 
inate the charge transport. Very low temperatures and 
clean samples are needed to expose it in experiments. 


Graphene stands out, among known materials, as 
the ultimate hydrodynamic material. Ultra-high-quality 
samples of graphene deposited on hexagonal Boron Ni¬ 
tride, which can be produced by standard techniques, are 
indeed virtually free of long- and short-range impurities 
and of charge inhomogeneities. 13 This, combined with the 
linear band dispersion and the lattice stiffness, 3 ' yields a 
fairly large temperature window in which the hydrody¬ 
namic regime can be realized. Indeed, due to the high lat¬ 
tice stiffness, scattering with in-plane acoustic phonons 
becomes relevant only at temperatures of the order of the 
room temperature 38 . 

Hydrodynamics is usually defined for Galilean- 
invariant systems 36 , in which the conservation of momen¬ 
tum leads to a diverging charge conductivity even in the 
presence of electron-electron interactions. However, at 
low energies graphene is well-described by a Dirac Hami- 
tonian, which is not Galilean invariant. This reflects 
in the fact that a homogenous current is not propor¬ 
tional to the total momentum, and can in principle be 
relaxed by electron-electron interactions, thus invalidat¬ 
ing the whole hydrodynamic picture. Our calculations 
show that, while the spin conductivity is finite and lim¬ 
ited by electron-electron interactions, the homogeneous 
charge conductivity is infinite as for an interacting 2DEG. 
The latter result stems from the fact that, at low temper¬ 
ature, doped graphene behaves as an effectively Galilean- 
invariant system, and paves the way to the description of 
charge transport in graphene in terms of Navier-Stokes 
equations. The relaxation of inhomogeneous current dis¬ 
tributions is controlled by the viscosities of the electron 
liquid in graphene, which will be given in a forthcoming 
publication. 

As the doping level decreases and graphene approaches 
the undoped regime, the momentum-velocity relation can 
no longer be linearized. Thus, not only the Drude weight 
gets renormalized, but also the charge transport time be¬ 
comes finite because of electron-electron interactions 39 . 
This in turn implies that, in this regime, the evolution of 
the charge current is described by a “generalized” Navier- 
Stokes equation which contains a relaxation term. While 
this Manuscript was being prepared, we became aware 
that a hydrodynamic theory of undoped graphene was 
worked out in Ref. 40. To describe the transport in 
this system in the collision-dominated limit, it is neces¬ 
sary to consider simultaneously the equations of motion 
of three macroscopic currents, i.e. the charge, energy 
and “quasiparticle-inrbalance” currents 40 . Note that the 
energy current, being proportional to the momentum in 
systems with a linear band dispersion, is always strictly 
conserved by electron-electron interactions and follows 
a Navier-Stokes equation. Conversely, the equation of 
motion of the charge current contains a finite relaxation 
rate due to electron-electron interactions. In the doped 
regime the three currents coincide 40 and the relaxation 
rate of the charge current due to electron-electron inter¬ 
actions vanishes. Thus, the transport can be described 
by taking into account only the charge current, whose 
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evolution is now described by the same Navier-Stokes 
equation satisfied by the energy current. 
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Appendix A: The calculation of the quasiparticle 
lifetime 


Before deriving the quasiparticle lifetime, let us recall 
Eqs. (11)-(13) which define the GW self-energy. They 
read 

T,^\k,ie n ) = - k B T y^ W(k' - k, ie n > - ie n ) 
k' ,A' £ n ' 

x G^{k',ie nl )V xx ,{k,k')V x , x {k',k) , 

(Al) 


and 


W(q,ifl m ) = 


1 (*7A^m) 


(A2) 


and finally 

Xnn(q,iu m ) = N v k B T Y E G[ a ,?(q\ie n ) 

q' ,E n ,cr' X " ,/J>" 

x Gy, 1 (q 1 + q, ie n + iu] m ) 

x / D\"i J ."{q',q' + q)^"A"(q / + q,q') ■ 

(A3) 


We first consider Eq. (Al), and we analytically con¬ 
tinue it to real frequencies. We thus define fs(i£ n ' — 
i£ n ,i£ n i) such that 


E {rT \k,i£ n ) = -k B T Y fc( i£ n' ~ i£ n ,i£ n ') 

E n t 

= j> ^n F ( z )fc{z - i£ n , z) . (A4) 

The contour C in the complex plane encircles all the poles 
of the Fermi function n F (z) = [exp {j3z) + l] , and 
leaves outside the branch cuts of fs(z — i£ n , z ), which are 
parallel to the real axis and pass through z = 0, i£ n . De¬ 
forming the contour of integration to encircle the branch 
cuts of fs(z — i£ n , z), and taking the limit i£ n —> £ + iq 
we obtain the retarded self-energy 14 

/ °° Hr 1 f 

— {[n F (0 +n B (£' -e)] 

X [Ms'- - e, e' + ) - /s(e' + - e, e' + )] 

+ n F (e') [/ E (e'_ - £, e'_) - /s(e' + - e, £+)] | . 

(A5) 


Here e± = e±iq. The term in the last line of Eq. (A5) is 
purely real, since it is multiplied by the imaginary unit. 
We thus get 

-i /»oo /*oo 

- = ~ 2 d£ 

1 ee J — oo J —oo 

x Y 3mG ( ' 7) (fc',e , + )3mW(fe-fc',e , + -e) 

x V xx ,{k,k')V x , x {k',k) . (A6) 

Here we understand that \k\ = k F . The imaginary part 
of the screened interaction reads 


de' dn F (£) 

7T ds 


[n F (£') + n B (s' - e)] 


%mW(q,u} + ) = \W(q,u) + )\$smXnn{q,u+) • (A7) 

The analytical continuation of x„ n (q, iw TO ) defined in 
Eq. (A3) can be performed in analogy with that of the 
self-energy. We thus define 


Xnn (*3S iu m ) = k B T Y^ fx(i £ n + i-Um,i£n) 

E n 

f dz 

= - j> ^Yi n F (z)fx(z + iu m ,z) .(A8) 

Also the contour C encircles only the poles of the 
Fermi function and leaves outside the branch cuts of 
fx( z + iui n ,z), which are parallel to the real axis and 
pass through z = 0, — ico m . Deforming the contour of in¬ 
tegration to encircle the branch cuts of f x (z + iui n ,z), 
and taking the limit ioj m —> to + ir\ we get 


Xnn(<?,W+) 


= ~J {[n F (e" + w)-n F (e")] 

x [/(e" + w, e"_) — /(e" + w, e "_)] 

+ n F (£") [/(e+ + u, e") — f(e'L + w, e "_)] | . 

(A9) 


Again the last term gives no contribution to the imagi¬ 
nary part, which becomes 

r ApH 

Smx„ n (g, uj + ) =-N v Y J — [™F( e ") ^ n F {e" + w)] 

<?V', 

x" ,/j," 

x SsmGy, - 1 (q 7 , s" )QmG^) ^ (q' + q, e" + w) 
x / D X "v,"(q',q'+ q)'Dn"y , (q' + q,q') ■ 

(A10) 


We put Eq. (A10) into Eq. (A7) and then back into 
Eq. (A6), and we get 



ds” <9n F (e) 
7r ds 


x [n F (e') + 71 b (e 7 — e)] [?i f (£ 77 + e) — n F (£ 77 + e 7 )] 
x | W{k - k^E^mG^ik^E^mG^iq' - fe,e" +e) 
x %mG%\q' - k',e" + e' + )V xx ,{ k, k')V x , x {k', k) 
x 2?A'V'(q 7 - k,q' - fc 7 )X> M " A "(q 7 - k’,q' - k) . (All) 
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We now use the fact that 

■A/ - = [ n F(e') + n B (s' - g)] 

x \n-F(e" + e) — ^(e" + e')\ 

= de" [ np(£ + £ ) _7lF ( £ )] 
x ^(e' + e") — np(e / )] 

„ 2 9nB(e") dn F {e) dnF(e') 
de" de de' 

In evaluating an integral of the form 

r °° ,dn B (£"),,2,///s 


X = 


de" 


cte" 


-e" 2 f(e") 


(A12) 


(A13) 


where /(e") is some smooth function of its argu¬ 
ment, we exploit the fact that the weighting function 
e" 2 dn B (e")/de" is strongly peaked at e" = 0 and its 
width scales with k^T 2 /e p. This does not mean, how¬ 
ever, that one can simply replace f{e") by /(0). Such 
a crude approximation would introduce a spurious diver¬ 
gence in the quasiparticle decay rate, because it spoils the 
subtle cancellation between two infinities which occur (i) 
in the phase space of the collinear scattering 41,42 and (ii) 
in the screening of e-e interactions. Both divergences are 
connected to the linear-in-momentum energy dispersion 
of massless Dirac fermions. The cancellation occurs as 
long as the argument of the function f(e) is finite. To 
take this into account we approximate 



de" 


driB^e ) ,,2 

de " A 1 


where e can be estimated as 


2t r 2 (fc B X) 2 
3 


/(e) + 0(T 4 ), 


(A14) 


£ 


1 / 3 f + °° ,, dn B (e) 

2 V 2t r2(fc B r) 2 /-oo de 


Te F . (A15) 


Here we have defined T = (k B T/eF and £ = ir/\/h. 
The factor — 3/[27r 2 (fcBT) 2 ] normalizes the weight of the 
function e 2 dnF,{e)/de to one. We have thus taken e to 
be half of the variance of the distribution e 2 dn B {e)/de. 
Eq. (A14) shows the crucial approximation that distin¬ 
guishes our results for the quasiparticle lifetime from 
those of Ref. 30. There the authors, although starting 
from the same GW approximation for the self-energy and 
deriving an expression equivalent to Eq. (All), did not 
approximate the final result according to Eq. (A14). The 
latter allows us to reduce the number of numerical inte¬ 
grations to be performed. With this approximation we 
finally get the quasiparticle lifetime at the Fermi surface: 

— = ~N v (k B T) 2 Y \W{k-k',e + )\ 

Tee k'.q'.a' 

x - fe,0 + ) 

x V - £+)£>++(&, k')V ++ {k', k) 

x X> ++ (q' — k,q' — k')'D ++ (q' — k ', q' — k) . 

(A16) 


Shifting k’ —> k — q" and q' —» k — k" we get 


— =~^N v (k B T) 2 Y IW,£+)I 

7”ee 3 inn, 

k" ,q" ,cr 


x $smG+\k — q" ,e + )^smG^ \k" — q",e + ) 

x SmG^tfc", 0 + )£> ++ (fc, k - q")V ++ (k - q", k) 
x V ++ (k",k" -q")V ++ (k" -q",k") , 

(A17) 

which can be recasted into the following Fermi-golden- 
rule form 

1 - 


qp 


3m[G A,7) (fc — q, —e) 


1 + COS (ifik-q - <Pk) 


(A18) 


Note that in these equation we do not sum over the spin 
index a. Since the system is spin-unpolarized the life¬ 
times at the Fermi surface of the two spin populations 
coincide. 


Appendix B: Analytical continuation of the response 
function 

To analytically continue ico m —! u, + iq in Eq. (9) we 
define the function f(ie, ie + iui) such that 

X-h V)(g>*w m ) = k B T Y f(ie n , ie n +iui m ) 

Ja Jp ' 

E n 

= ~ ^n F (z)f(z,z + iu m ) , 

(Bl) 

where we suppress for brevity all its momentum depen¬ 
dence. The contour of integration in Eq. (Bl) is chosen 
in such a way as to encircle all the poles of n F {z) and 
to exclude all the branch cuts of f(z,z + iui m ), which 
occur for ^sm(z) = 0 and £ sm(z + iui m ) = 0. We trans¬ 
form the contour-integration in an integration over the 
branch cuts and we then perform the analytical continu¬ 
ation ioj m -A u! + iq. We get 

w ) = - J { [«f(£ + w) - n F (e)] 

X [/(£_,£+ +W) -/(£_,£_ +W)] 

+ n F (e) [/(£+,£+ +w) — /(£_,£_ +w)] | . (B2) 

Here e± = e ± iq. After the analytical continua¬ 
tion, G^\k,e + ) = G^' a ' 1 (k,e) and G^(fc,£_) = 
G^ ,cr \k,e). Here G^’ a \k,e) [G^ R,<T ^(fe,£)] represents 
the advanced [retarded] Green’s function. 
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Note that the square brackets in the last line of 
Eq. (B2) contain a purely imaginary quantity, which 
(being divided by the imaginary unit) gives a purely 
real contribution to x -O) -o') (Qi u )- Note also that 

Ja Jp 

f {£-,£+ A co) contains the product of a retarded and an 
advanced Green’s function, whereas in /(e_,e_ Aw) and 
/(e+, £ + Aw) both Green’s functions are either advanced 
or retarded. The last two functions [/(£_,£_ + w) and 
/(£+, £ + + w)] have all the poles on the same side of the 
complex plane. Note however that, as usual 1 * 2 , we can 
exploit this property only performing the integral over 
the band energies. We thus get that /(£_,£+ Aw) gives 
the dominant contribution when £pr ee 1- Since we are 
interested in the Fermi liquid regime, in what follows we 
will retain only this term. Eq. (B2) thus becomes 


an integration over the branch cuts of g(z, z A ico m , z — 
i£ n ). We then perform the analytic continuations with 
the prescription iio m —> w+, i£ n -A £_, i£ n A ico m -A 
£+ A co. After some lengthy but straightforward algebra 
we get 

A (1 ’ 2) (fe + ,£+ + W, fc-,£_) = J ^{n F (e') 

X [g(e' + , £+ + w, e' + - e) - g (£'_, e' + A w, £+ - e)] 

A np^ 1 ) \g{£'_ — W, e' + , e'_ — £ — w) 

- g{£'_ - w, e '_,£'_-£- w)] 

— ns(£ / ) \g{^'- A £> £+ + £ + w, e' + ) 

~ g( £ - +£,£+ + £ + w ,£-)] | ■ (C2) 


/ de 

— [n F (e Aw) -n F (e-)] 

K,A,A ^ 

X Gi A ’ CT) (fe_, £) A&f (fc_, (fc + , £ + w) 

x A yxM k +> £ ++ 0J , k -,£-) , ( B3 ) 

which coincides with Eq. (23). 


Appendix C: Analytical continuation of the 
Bethe-Salpeter equation 

In this section we guide the reader through the long 
and complicated calculation of the vertex correction to 
the thermal conductivity. We analytically continue the 
three contributions AvX^(fe+,t£ n + ico m , fc_, i£ n ) (i = 
1,..., 3), defined in Eqs. (15)-(19) to real frequencies. In 
the Fermi-liquid regime we consider only the dominant 
contribution to the dressed vertex, to be used in com¬ 
bination with the product of the retarded and advanced 
Green’s functions that appears in Eq. (B3). From this 

we see that the analytic continuation of ^■y’x'g ^ ( k + > ® £ n + 
ico m , k_,i£ n ) is done with the prescriptions ico m -A w+, 
is n —y £~, iSn + ^ £-\- H - 


1. Analytical continuation of Eq. (15) 

We define the function g(i£ n *,i£ n t A ico m ,i£ n ' — i£ n ) 
such that 


We now shift £' —> £ l + co in the third and fourth lines of 
Eq. (C2). We note that we can safely take the limit w —>• 0 
in nF(e , + w). Note also that g(e' + ,e' + A co,s' + — e) and 
g{e'_ — co,£'_,£'_ — £ — co) have the poles on the same side 
of the complex plane, and therefore can be neglected 32 
in the limit of £FT ee 1. Shifting e' —> d + £ in the last 
two lines of Eq. (C2) we readily obtain 


A (ll2) (fc+,£ + Aw,fc_,£_) = ^ Y / ^ 

k',cr" 

X (k\ k, £'_ -£)- (k>, k, e' + - £) 

x [n F (e') + n B (e' - s)]G^ ,(T \k' + ,e' + u) 
xG< A/,) (fc',£') A ( ;'${k' + ,£' ++U ,k'_,£'_) . (C3) 


It remains to determine W^y’^7 \k’,k,£± — e). 
Eq. (C3) implies that we have to analytically continue 
the functions W* 1,2 ' 1 with the prescription ie n > -A e'_ and 

Z£ n ' A ZCO m —^ £(|_ A W. 


a. The analytical continuation of Eq. (17) 

We now perform the analytical continuation of Eq. (17) 
with the prescription *w m -A w+, i£ n -X £_, i£ n A ico m —t 
£_j_ A co, i£ n i -A £ r _, and ie n i A ico m e + A co. As 
shown in Eq. (C3), we need to calculate the function 

W xx^'\ k ^ k 'd- “ £ ) “ W xx^'\ k ^ k ^ e '+ - £ )’ which 
reads 


A*- ^ (^+5 ^£n “I - 5 ^n) 

= ^ ^ £l('i'£n f 5 “I - ^ £n r ^n) 

£ n r 

/ dz 

n F {z)g(z, z A iu m , z - i£ n ) ■ (Cl) 


Here and in what follows we suppress for brevity all the 
spin, band and spatial indices of the dressed vertex. As 
usual, we transform the sum over the poles of n-p(z) in 


W&ffi (*b - *) - W^\k, k’, £ ' + - £) = 

x [W(k — k', e'_ — e) — W(k — k', e' + — e)] 

x fc + )I? M/ \(/c_, fc_) 

= — 2i5 aa "^sm[W (k — fc', £' + — e)] 

x 'Dx’n' {k+, k + )D il x{k _, fc_) 

= ~2i6 aa „\W{k - k',e' - £)| 2 3m[xi°2(fc - k',e' + - e)] 
x P AV (fe+,feV)^A(feL,fe_) ■ (C4) 
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The imaginary part of the density-density response func¬ 
tion was given in Eq. (A10). After some straightforward 
manipulation, we get 


dui' 


i[xnn(k- k',e' + -e)] =2N V ^ E / 


q' ,a"' A",/^" 


x + e') — uf(uj' + £)] 

xSm \q — k,uj' + e) 

\q' — fc 7 , uj' + e') 


x^sm 


x£>a>"( q' — k,q' — k')V^ X " (q' - k'. q' - k) (C5) 


Putting Eq. (C5) into Eq. (C4) we finally find 


W^\k,k',e'_ - e) - W^\k,k',e' + - e) 

= 4N v S a(J "\W(k - k! , £ 7 - e)\ 2 E E 
x J [ n F (w 7 + e') - n F (w 7 + e)] 


x3m 


,(R,o-" 

A" 

V 

fi" 

V 


x2?av' (^-H )*Dn\(k— i k — ) 

x£>a>"( q' — k,q' — k')V^ X " (<?' - q' - k) ■ 

(C6) 


b. T7ie analytical continuation of Eq. (18) 


We now turn to the analytical continuation of Eq. (18) 
with the prescription iuj m —i ui+ , ie n —> e_, IE™ + iw m —>• 
£+ + ui, ie n > —> £ 7 _, and ie n i + ico m —► £ 7 + + w. This time 
we define 


fc A £ n' - i£ n) = j> ^«b(z) 

X tU 2 (2, 2 - *£„ + - 2, ie n > + iu) m - 2:) . 

(C7) 


Integrating over the branch cuts of W 2 (z,z — iuj m ,i£ n + 
iui m — z,i£ n > + iu} m — z), and performing the analytical 


continuations as stated before, we get 

] (4 fc ’ 4 - £ ) = J 7^1 {^b(w') 

x [«; 2 (u/|_,u/_ — uj, £ + uj — uj'_,e' + w — ui'_) 

— W2 (u/_ , U>'_ — W, £ + UJ — Uj '_, e! + UJ — Uj'_ )] 

+ n B (u/) [zy 2 (w+ + <*;,<*;+, £ — o/|_, £ 7 - (J + ) 

- W 2 (oj' + + UJ, uj'_,£ - Uj' + ,e' — u/)_)] 

— n F (ti/) [tt> 2 (u/,_ + £ + W, + £, — u/|_, £ 7 — £ — w(p) 

- w 2 (w+ + £ + W, w'_ + £, -Uj'_,e' - £ - w(p)] 

- nF(w')[w2(u' + + £' + W, + £ 7 , £ — £ 7 - — W+) 

— W 2 (w(_ + £ 7 + CC, + £ 7 , £ — £ 7 — Ca(j_, — Uj'_ )] . 

(C8) 

Note that the terms on the r.h.s. of Eq. (C8) proportional 

(c\ t/\ 

to nB(w') are identical in both W^?° (fc 7 , k , e'± — e) and 
thus vanish when the difference is taken. We thus neglect 
them in what follows. Eq. (C8) reduces to 

fc -4 - £ ) = ^ J ^“{ nF 4' - £ - w ) 

X [w 2 (w(|_, U}'_ — UJ, £ + UJ — o/|_, £ 7 + UI — w(p) 

— W2(oj ' + , U)'_ — UJ, £ + UJ — <*/_, e' + UJ ~ 4)] 

+ uf(uj' — e' — to) 

x [rc 2 (cc_j_ A 7 oj_ — uj i £ -\- uj — w_j_, £ 4 - w — ) 

— W 2 (w(|_, U)'_ — UJ,£ + UJ — UJ± ,e'+uj — u>'_ )] | ■ 

(C9) 


Finally, 


IE 


(2,era 




i > (£'_~e)-W^ u ,\e',-£) 


d&J 

27TZ 


X 

X 


X 

X 

X 

X 

X 


[np(c</ — e — cu) — — s’ — a;)] 

[w2(w'+, 0j'_ — UJ, £ + UJ — 0j' + ,£ f + Id — LJ_ 
W2 , (jo'_ — 0J,£ + UJ — , £ f + UJ — Lj'_) 

W2(id , + ,UJ , _ — i0,£ + UJ — a/j_, £ r + UJ — UJ+) 
W2 (lo'+ , to'_ — (J, £ + CJ — (jj'_,£ r -\- OJ — LU+) 

ANv E / ^(q>;)W(q>_- 


) 


w) 


[?rF(w / — £ — aj) — nF(w 7 — £ 7 — w)] 
3m[G^:’ CT ^(fc + — q 7 , £ + w — a/)] 

3m [G^’ 17 ^ (&+ - q 7 , £ 7 + uj - w 7 )] 

2?A'A"(fc+,fc+ - q')V X "\{k + - q',k _) 

(fc'_, fc; - q')^"/,' (K — q'i k +) ■ (cio) 
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We can now take the limit v F q «w« e f , and we get manipulations Eq. (C13) becomes 


W- 


(2,era") , 


XX' 


/ } ( £ - - } (4 - £ ) = 4^Vv E 


X / W ^| 2 [ nF ( w/ “ £ ) - ™f( w ' - e')] 

x3m[G^ ,<T ^(fe — q',e — a/)] 

x3m[G^ ,<T — q', e'— a/)] 

xVyx"(k,k - q')Dx"x{k - q',fc) 
xD w » ( fc/ ’ k ' ~ ( k' - q\ k’) . 

(Cll) 


A( 3 )(fc + ,£ + + W ,fc_,£_) = E E / |^ 

k',cr" /x,/x' 

X W^,\k,k',£'_ +£) - w^\k,k',e' + +£) 
x [n F ( e 7 ) +«b(£' + e)]G^’ ct ) (fcV> £/ + w ) 
xG^"Hfc'_,£OAE?(fc;,£ ++W) fc'_, e '_) . (C14) 


It only remains to determine ^ (k' — fc,£± — e), 

defined in Eq. (19). Eq. (C14) implies that we have to 
analytically continue the functions HA 3 -* for ie n i —> e'_ 
and is n ' + ito m —> e' + + to. 


2. Analytical continuation of Eq. (16) 


We define h(ie n ',ie n r + ico m , is n > + i£ n +ico m ) such that 


The analytical continuation of Eq. (19) 


A ( 3 \k + > iErn 4“ i k -, t£n) 

= kj$T A ) h(is n <, is n ' T iuj rn , is n ' -f- 7£ n T tojm) 

e n t 


clz 


27 ri 


-n F (z)h(z, z + ico m , z + ie n + ito m ) . 


(C12) 


We now turn to the analytical continuation of Eq. (19) 
with the prescription iu} m —> w+, ie n —> £_, ie n + ibj m 
£+ + to, ie n t —> e'_, and ie n ' + ito m —)•£(,_+ to. We define 


W xx qE( fc, ’ fc ’* £n ' +i£n + iu} m) = j> ^An B (z) 
x w 3 (z,z - iio m ,i£ n + z,i£ n '+ ito m - z) . (C15) 


Here and in what follows we suppress for brevity all the 
spin, band and spatial indices. To perform the analytical 
continuation we first transform the sum over the poles of 
n F (z) in an integration over the branch cuts of h(z,z + 
ito m , z + i£ n + *w m )- We then analytically continue the 
result, according to the prescription ioj m —> to+, is n —> 
£_, i£ n + ito m —)•£++ a;. After some lengthy algebra we 
get 

A (3) (fe+,£+ +w,fc_,£_) = J ^{n F (£') 

X [h{£' + , e(|_ + 10, e' + + £ + to) 

— /i(e(_ , e\ + W, e(j_ + £ + w)] 

+ n F (£ / ) [/i(e(_ — 00 , e' + , £ ; _ + e) 

— h(e'_ — w, e'_, e'_ + e)] 

— 7Zb(£ / ) [M e - — £ — W, £+ — £, £+) 

— }i(e'_ — £ — W, £+ — £,£(_)] | . (C13) 

We now shift e' —► e' + w in the third and fourth lines of 
Eq. (C13), and we take the limit to —> 0 in ti f (e' +to). We 
note that /i(e+, e' + +io, £++£+<*;) and h{E > _ — to, e'_,e'_-\-£) 
have the poles on the same side of the complex plane, and 
can be neglected in the limit £ F r ee 1. We then shift 
e' —X e' + £ + <j 0 in the last two lines of Eq. (C13), and 
we take the limit to —► 0 in ^(e 7 + £ + u>). After these 


Integrating over the branch cuts of w 3 (z,z — ioo m ,i£ n + 
z,i£ n > + iuJm ~ z), and performing the analytical contin¬ 
uations according to the prescriptions stated before, we 
get 

w \\W } ( fe '> fc, e± + £ + w) = J {n B (w') 

X \w 3 ((jo' + ,(jo'_ — UJ,Co'_ + £, e' + 00 — U)'_) 

— W 3 ((jo'_, Uj'_ — UO,Uo'_ +£,£'+ UO — w(_)] 

+ 71b [ w 3 ( w + T , w!)_ + £ + UO, e' — w(|_) 

— w 3 (uj' + + UJ,Oo'_, uj' + + E + 00 , e’ - <*/,_)] 

— n F (a; , ) 

X [w 3 (oo' + — £,io'_ — £ — 00, lo' + , e' + E + 00 — o/p) 

— W 3 (lo' + — E,<J_ — £ — UJ,0o'_,e' + £ +07 — w(p)] 

— n F (o/) 

X [ui 3 (w(|_ + £ 7 + 00, u/_ + E 7 , io '_|_ + £ + e' + 00, —w(j_) 

— W 3 (lo'_i_ + £ 7 + 00, U}'_ + E 7 , + £ + £ 7 + CO, —C*/_)] . 

(C16) 

Note that the terms on the r.h.s. of Eq. (C16) propor¬ 
tional to 71b {to') are identical in both ,\k' , k,E'± — 

e) and thus vanish when the difference is taken. We thus 
neglect these terms in what follows. Eq. (C16) reduces 
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to 


w \\ ^V( fe, ’ fe > e ±+ e+w ) = - J ^{mw'+e) 

x [ 103 ( 0 ;+, J_ — to, u/|_ + £, e’ + w — uj' t ) 

— W 3 (u)' + , w(_ — CO, U}'_ + £, s' + CO — C0+)] 

+ 71f(w / — e' — Co) [l 03 (c0+, CO*) — CO, C0+ + £, £ 7 + CO — C0+) 

- W 3 (uj' + ,U}'_ - CO,CO+ + £,s' + CO - Co'_)] | . 

(C17) 

Finally, substituting Eq. (19) into Eq. (C17) we get 


+*+")- + £ +“) = / D 

x [77+(cu 7 + e) — nF(u/ — e’ — co)] 
x [to 3 (co+,co(_ — co, a/, + £, £ 7 + co — 

— w 3 (u' + , lu'_, co(_ + £, e ' + co — co(_) 

- w 3 (uj' + , co 7 _ — u>,u}' + + £, £ 7 + co - co+) 

+ u> 3 (co+, co 7 _ — co,u/_ + £, £ 7 + co - co+)] 

= -4 N v EE / 

q' A ",y." J 

x [?7f(co 7 + e ) — ?7f(co 7 — s ' — co)] 
x 3m[G^ ,1T *(fc_+q',£ + w')] 
x 3m[G^’ ff ^(fc+— q',£ 7 + co — co 7 )] 
x X> AA //(fc_,fc_ + q')V X "\'(k^ +q',k+) 
x 2? MM »(fe , _,fcV-q , )7 ? M'v( fc +-9',fcV) ■ (C18) 


Taking the limit opq <u<£f, Eq- (C18) becomes 


Ma’vm'V- + £ ) - W K'V+ + £ ) = -4JV v X! 

q',a" 

x E / T“l M/ (9'> w, )| 2 hF(w / + e) -n F {u' - £ 7 )] 

\",n" •> m 

x 9m [G^ ,<t) (fc + q 7 , £ + co 7 )] 
x9m[G^’ CT \k 7 — q',£ 7 — co 7 )] 
x 77 AA " (fc, k + q')V X "\’ (k + q 1 , k) 


xV^n(k',k' - tftV^^k! - q',fc 7 ) . 


(C19) 


Appendix D: The derivation of the Bethe-Salpeter 
equation in the charge and spin channel 


Using Eq. (29) in Eqs. (25)-(27) we get 

8iN v 


A^\k,e ++ .,k,e-) = - 


co + i/r e , 


E E 

k' ,q' ,a" n,fi" 


x [nF(w 7 + £ 7 ) — 71f(co 7 + £)] \W(k — k',e' — e)| 2 
x 3m[G^ ,<T \q' — fc,co 7 +£)]3m[G[ l R ’^(fc 7 ,£ 7 )] 
x Sm[Gfi a "\q' - fc 7 , J + e')]V x ^k, k') 
x 'D fjX (k', k)T> X n (q 7 — fc, q 7 — k') 
x V^x"(q' - k',q' - k)A^J(k', e' + + co, fc 7 , £ 7 _) , 

(Dl) 

and 


A%j\k,e ++U ,k,e-) = - 


8iN v 


co + t/r e , 


E E 

k' ,q' ,cr" /i,, X" ,/i" 


X + »=(£'-£)] 
x [?7f(co 7 — e) — uf(w 7 — £ 7 )] |IU(q 7 , co 7 )| 2 
x Sjm[G|^ ,cr '(fe 7 — q',£ 7 — w')]3m[G| 1 R ’' 7 ^(fc 7 ,£ 7 )] 

x 3?n[G^’ CT) (fc - q 7 ,£ - w 7 )]2? AA "(fe, k - q 7 ) 
x V X n X (k - q’,k)V^n{k’, k 1 - q , )'D t _ l ” l _ l (k' - q 7 ,fc 7 ) 


Eu./ 3 ) ( fc, ^ £ + +W,fc 7 ,£ 7 _) , 


W»,/3 

and finally 


(D2) 


A AaX } ( fe ’ e + + W,fc,£_) = 


8iN v 


+ i/ 7 


■ E E 

k' ,q' ,cr" fi,X" ,fi" 


x / £/£[’*(£) + »Biy + r)] 
x [n F (w 7 + s) - n F (w 7 - e 7 )] |lU(q 7 , w 7 )| 2 
x 3m[G^’°^(fc + q 7 ,w 7 + e)] SrofGjf'" (fc 7 , £ 7 )] 

x 3m[G^ i ;’ <T \k' — q',£ 7 — u/)]2? AA "(fc, k + q') 
x V X n X (k + q 7 , k)T> (fc 7 , k! - q')T> M - v (fc 7 - q 7 ,fc 7 ) 
x - (D3) 

In these equation the limit u> —> 0 is understood. 
Eqs. (D1)-(D3) should be plugged into Eq. (24) and then 
back into Eq. (B3) taken in the limit oj —> 0. The latter 
is given in Eq. (28), and can be further approximated as 


, _ n v 2 iN v 

XjMAr'nQ — 0,w) —>■ — uj ——— 

Jo, Jp uj 4 - %t 


?/ 


de dn-p{e) 


2m de 

x Sm[Gf | R ’ o ’ ) (fc,0)]A^ ) a (fc,fc)A^J(fc,£ + +w,fc,£_) . 

(D4) 
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To obtain this expression we used Eq. (29), together 
with the fact that the derivative of the Fermi function 
is peaked in e = 0, and that 9m[G A R, ' T (fc, 0)] constrains 
|fc| = kp and A = +. We now approximate Eqs. (Dl)- 
(D3) by noting that the combination of Fermi and Bose 
functions constrains £ ~ e' ~ uJ ~ 0. We thus get 


A (1><™F) 

++,/3 


(fe,e+ + w, fe,e_) 


8 iN v 
uj + i/r ee 


E 

k' ,q' ,a" 



x [n F (u/ + e 7 ) - n F (w' + e)] |W(q 7 , 0)| 2 
x - g',0)]9m[G^ ,,T ^(fe 7 ,0)] 

x 3m[<E R,,T " V - k 1 , 0 )]V++(k, k - q')V ++ {k 
x V ++ (k',k' -q')V ++ (k' -q',k') 
x - q',w + ,k - q 7 ,0 _ ) , 


q',fc) 


(D5) 


together 

3 


E A ++E fc E,M ) = ~; 


8iJV v 


+ i/7 


E iw,o)i- 




x EE Me'l+nele'-e)] 

x [tif(u/ + £ 7 ) — 71f(w' + e)]Sm[G^ ,,T ^(fc',0)] 

x $sm[G^’ a \k — q',0)]9m[G+ l ’ <T \q' — fc',0)] 
x V ++ (k,k - q')V ++ (k - q',k) 
x 2? ++ (fc 7 ,fc 7 -q 7 )2) ++ (fc 7 -q 7 ,fc 7 ) 

x A++i( fe - w + , fe - q', 0”) + A^^^fe', w + ,fc',(T) 

- A2.' <7 ' ) (fc , -g , ,o;+ fc'-g'.O-)! . (D8) 


Plugging Eq. (D8) back into Eq. (24) we immediately get 
the self-consistent Bethe-Salpeter Eq. (30). 


and 


Appendix E: The spin transport time 


» (2,crcr') 

++,/3 


(k, £ + + OJ, k, £_ 


) 


8 iN v 
u + i/r ee 


E 

•/ n' n 



x [np{J + e') - n F (J +e)] |W(q 7 ,0)| 2 
x 9m[G^| R,CT " ) (A: 7 - q 7 ,0)] 9m [Gf 0 (fc',0)] 
x $sm[G^’ a \k — q 7 ,0)]P ++ (fc, k — q') 

X V ++ {k - q\ k)V++(k', k' - q')V ++ {k' - q 7 , k') 
x A^^fe',^,fc 7 ,°-) , (D6) 


and finally 


a (3,<tct') 

A ++,/3 


(k,e + + w,fc,e_) 


8 iN v 
w + i/r ee 


E 

k',q',a" 



x [n F (w 7 +e') - np(J + e)] |W(q 7 ,0)| 2 
x 9m[G^' ff, (fc - g',0)]9m[G ( J !,a \k' — q', 0)] 

x 9m[G^ R ’ <T " ) (fc 7 ,0)] Z?_|_|_ (fe, k — q 7 ) 

x P++(fe - q 7 , fc)T» ++ (fc 7 - q 7 , k')V ++ (k', k! - q 7 ) 
x A^ + ^\k'- q',uj + ,k' -q',0 - ) . (D7) 


We shifted k! —> k — q' and q 7 —> k — fe 7 in Eq. (D5), 
u/ -A -w 7 in Eq. (D6), and e 7 —> —s', k’ — » k’ — q 7 
and q 7 —► —q 7 in Eq. (D7). Moreover, we used that 
the imaginary parts of the Green’s functions constrain 
their momentum argument to the Fermi surface and the 
band indices to be all equal to Putting everything 


The spin transport time defined in Eq. (37), which we 
recall here for completeness, reads 

4) = ~lN v (k B T) 2 J2 |^(q',0)| 2 9m[Gf ) (fe —q 7 ,0)] 

^tr k',q' 

x 9?n[GV R) (fc 7 ,0)]9m[G^ R) (q 7 - fe 7 ,0)] 
x V ++ (k,k - q')V ++ (k - q\k) 
x 2? ++ (fc 7 ,fc 7 -q 7 )2) ++ (fc 7 -q 7 ,fc 7 ) 
x [1 - cos(tp k - ipk-q>)\ ■ (El) 

In the low-temperature limit the three Green’s functions 
on the right-hand side of Eq. (El) can be approximated 
by (5-functions, which constrain 


cos (ip k - <p q >) = 
cos (ipu - <p q >) = 


E 

2/cp 

E 

2/cp 


(E2) 


and thus 1 — cos {ipk — Vk-q') = q' 2 / (2fcf). The Coulomb 
interaction forces us to disregard the solution of the three 
(5-functions with q' = 0. 

Note that Eq. (El) describes the simultaneous scatter¬ 
ing of a particle from the state k to the state k — q', both 
at the Fermi surface, and the creation of a particle-hole 
pair with total momentum q 7 . According to Fig. 3, the 
angles between k and k—q 1 , and k' and k' — q', are identi¬ 
cal. This implies that U + +[k, k — q') = V ++ (k\ k' — q 7 ). 
Moreover, 

*>++(*, k - q')V ++ (k - q 7 , fe) = 1 - (J-'j . (E3) 
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Putting everything together we get 

4) = -y^v(feT) 2 ^ |lP(q',0)| 2 3m[GV R) (fc-q',0)] 

^"tr k' ,q' 

x 3m[GV R) (fc , ,0)]3m[G^ R) (q' - fc',0)] 



Shifting w -> (p k > + ip k and ip q i —> ip q > + ip k , we im¬ 
mediately see that the three ^-functions are solved by 
ipq' = iarccos [q'/(2kp)\ and tp k ' = 0,2 <p q i. All these 
solutions give identical contributions. Summing all of 
them we finally get 


» 


N v {k B T) 2 r 2kF 
37ru|e F J o 


dq’ q'\W(q',0)\ 2 


i-(i- 

\ 2 kp 

2nN v a 2 ee {k B T ) 2 
3 hsp 

4(1 — 3_/V 2 a 2 e )arccoth(l + 27V v a ee )] . (E5) 


[3(2iV v a ee - 1) 


properties at the Fermi surface, we will consider in what 
follows the following one band model of graphene 

'U = Y^v F p i + ^^2 v ( r i - r j ) , (F2) 

* i¥^i 

where i and j label particles in conduction band. The 
interaction with states in the valence band provides the 
well-known logarithmic divergence 15 of the renormalized 
Fermi velocity Vp. We now imagine to perform the fol¬ 
lowing unitary transformation 


U = exp 



(F3) 


on the Hamiltonian of Eq. (F2). In Eq. (F3) we intro¬ 
duced the Pauli matrix t z which acts on the spin degree of 
freedom. Clearly the transformation U commutes with 
the interaction Hamiltonian, since it contains only the 
position operator r,;. To first order in q it induces the 
change in energy A E = j s ■ q. For the state |p,f), with 
the extra quasiparticle of momentum p and spin up we 
thus get 


Here we used that for 0 < q' < 2kp W(q', 0) reduces to 
the Thomas-Fermi interaction, with qtf = 2JV v o:ee&F as 
the screening wavevector. In the limit a ee —> 0 we get 


1 



AnN v ( k B T ) 2 
3 Hsp 


a ee ln(A v a ee ) . 


(E6) 


The logarithmic dependence on the coupling constant a ee 
is due to the Thomas-Fermi screening. 


A E = v s p ■ q . (F4) 

On the other hand, we may consider the variation of 
energy due to the shift generated by the unitary trans¬ 
formation of Eq. (F3) on the phenomenological energy 
functional of Landau’s theory of normal Fermi liquids 14 . 
Shifting the momentum of spin up (down) particles by q 
(-<?) we get 


Appendix F: The spin velocity 

We now show a brief derivation of the renormalization 
of the spin Drude weight of Eq. (41), based on Landau 
theory of normal Fermi liquids. The derivation closely 
follows that of Ref. 18. The question we answer in this 
appendix is: what is the spin current carried by a quasi¬ 
particle? We thus consider a state in which a quasipar¬ 
ticle with momentum p and spin up is added to the sys¬ 
tem. The adiabatically turned-on electron-electron inter¬ 
actions dress the quasiparticle and renormalize the spin 
current it carries. We consider the expectation value of 
the spin current operator j s = j- j- — on this state. In 
analogy with the charge current, we know that the spin 
current is proportional to the unit vector p. However, 
the constant of proportionality, namely the spin velocity 
v s , is to be determined. We define v s as 

(P,t |js|p,t) = v s P ■ (FI) 

where |p,t) denotes the full nrany-body state with the 
extra quasiparticle with momentum p and spin up. 

We now connect the phenomenological theory of Lan¬ 
dau to the microscopic model. Since we are interested in 


A E = vpp ■ q - ^2frt yP ' T Tq ■ V p >n 0 , T (p') 

P' 5 T 

= «f(1 + F?)p-q, (F5) 

where f P a. P ' T is the Landau interaction function 14 and 
n o ,t(p') is the equilibrium distribution function of quasi¬ 
particles. Note that at, if p and p' are both at the Fermi 
surface f pa .p'T is a function only of the difference between 
the angles of p and p’ , i.e. f pcr . p > T = -<P P >)- The 

Landau parameters F r (, with ^ = a, s, are defined as 14 

= 2^ / S [ /ttW ± /ti(0)] C ° S2(0) ’ (F6) 

F J 

where the plus (minus) sign in square brackets on the 
right-hand side of Eq. (F6) holds for t = s (£ = a). 

A direct comparison of Eq. (F4) with Eq. (F5) imme¬ 
diately gives 

v s = v£(l + F?). (F7) 

We obtain the renormalization of the spin Drude weight 
by replacing v s in lieu of the bare Fermi velocity vp in its 
non-interacting expression. 
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